Loschmidt echo in the Bose-Hubbard modeh turning back time in an optical lattice 
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I show how to perform a Loschmidt echo (time reversal) in the Bose-Hubbard model implemented 
with cold bosonic atoms in an optical lattice. The echo is obtained by applying a linear phase imprint 
on the lattice and a change in magnetic field to tune the boson-boson scattering length through a 
Feshbach resonance. I discuss how the echo can measure the fidelity of the quantum simulation, 
the intensity of an external potential (e.g. gravity), or the critical point of the superfluid- insulator 
quantum phase transition. 



PACS numbers; 03.75.Lm,39.25.+k,32.80.Qk 

More than two decades ago 0, Feynman envisioned 
quantum physics simulations being performed by control- 
lable quantum systems as no traditional computer could. 
For large complex systems, it is not always practical to 
estimate the accuracy of the simulation, given by the 
fidelity f{t) of the actual experimental evolution with re- 
spect to the desired one. Yet, in some cases f{t) can 
be measured directly. Let Tisim be the Hamiltonian to 
be simulated and Ti. = TLsim + V the laboratory Hamil- 
tonian, where V contains unknown terms present in the 
real experiment. To measure f{t), one has to first evolve 
a given initial state |'0o) with TL, perform an operation 
that changes the sign of H-sim, evolve for another time 
t with the new Hamiltonian Ti = —Tisim + |3> ^.nd 
then measure the probability of finding the system in the 
initial state. 
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Clearly, /(i)|y^o ^ 1- general, f{t) - also known as 
a Loschmidt echo (LE) - depends on V 3]. Therefore, 
to measure f{t) one only needs the ability to change the 
sign of the Hamiltonian. Notice that, in quantum me- 
chanics, this operation is equivalent to reversing time. 
Thus, the LE is the probability of return to the initial 
state by a forward-backwards evolution in presence of a 
perturbation. 

Recently, the Bose-Hubbard model (BHM) was simu- 
lated using cold bosonic atoms loaded in an optical lattice 
Many models can be simulated with this system 
However, the BHM has a strong appeal because of its 
simplicity and because it presents a superfluid-insulator 
transition - a paradigm of quantum phase transitions 
. fn this letter 1 propose an experimental procedure to 
change the sign of the BHM Hamiltonian and perform a 
Loschmidt echo. The time reversal operation consists of 
a phase imprinting in the lattice and a sign change of the 
boson-boson scattering length by varying the magnetic 
field near a Feshbach resonance - both techniques have 
been demonstrated 0,0,0- Also important, when the 
fidelity of the simulation is high, the LE can be a sensing 
tool: f will show how it can be used to measure e.g. the 
intensity of external potentials, or the critical point of 



the BHM quantum phase transition. 

A simple example of a LE is the Hahn or spin echo 
In a typical nuclear magnetic resonance (NMR) experi- 
ment, the polarization signal decays rapidly due to the 
inhomogeneity of the local magnetic field - spins precess 
at different Larmor frequencies. However, a tt rotation of 
the spins changes the effective sign of the spin-magnetic 
field interaction, refocusing the polarization as if time 
were reversed. The spin echo decay is used to measure 
the relaxation time T2, given by other terms in the full 
Hamiltonian (e.g. spin-spin interactions) that are not re- 
versed by the tt rotation. More complex LEs have been 
performed e.g. in solid state NMR [ll| (including many- 
body interaction terms), and in trapped cold atomic sys- 
tems pjj. In general, the LE is a measure of the stability 
of quantum evolution with respect to perturbations 
For instance, in a classically chaotic system there is a 
regime where the the LE decays with the Lyapunov ex- 
ponent of the classical Hamiltonian Also, in open 
systems the decay rate of the LE equals the decoherence 
rate Recently, it was shown that the decay of the 
LE can si gna l a quan tum phase transition in the environ- 
ment [ij, lla, lid ■ Foi' pure states, the LE is equal to 
the fidelity of a quantum computation fl? * 

The BHM Hamiltonian is 
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where are bosonic annihilation operators in the 
site of a discrete lattice, J is the hopping amplitude be- 
tween neighboring sites, and U is the interaction energy 
of bosons in the same site. The BHM undergoes a quan- 
tum phase transition when J/U is varied For J ^ U, 
the dominating hopping term favors delocalization of par- 
ticles: the ground state is a superfluid. In the opposite 
regime, J <^ U, the interaction energy is too strong and 
number fluctuations in a site are costly: the ground state 
(for integer density) is a gapped Mott insulator. The 
BHM is implemented using cold atoms in a periodic oscil- 
latory potential created with a standing wave laser light 
(the optical lattice) 0. Up to good approximation, the 
atoms occupy only the ground state of each well of the 
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lattice. The overlap between ground states of adjacent 
sites gives the tunneling amplitude J . Longer neighbor 
tunneling is suppressed exponentially. A contact inter- 
action potential with a range much smaller than the size 
of the wells is assumed. In the right limit, this trans- 
lates in the interaction term of Eq. Q , with U directly 
proportional to the s-wave boson-boson scattering length 
as- In all implementations of the BHM in optical lattices 
so far there is an additional magnetic trapping potential 
. To perform a LE it would also be necessary to change 
the sign of this potential (the system remains stable for 
a time equal to the forward in time evolution). However, 
an inhomogeneous magnetic field implies a position de- 
pendent interaction strength [/, departing from the sim- 
plest BHM. For simplicity, I will only consider the homo- 
geneous BHM, which will be realized experimentally in 
the near future 0,1131 ■ In addition, I will only consider 
systems with open boundary conditions. The experiment 
proposed here could be applied to other implementations 
of the BHM, e.g. with trapped ions or even in the 
Hubbard model with fermions. 

The change of sign of TIbh is done in separate steps 
for J and U . The change of U is achieved using a mar- 
velous experimental handle of cold atomic systems: the 
effective atom-atom scattering length can be tuned by 
varying an external magnetic field near a Feshbach res- 
onance 0. The external field shifts the energy of the 
free atoms with respect to a bound molecular state that 
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FIG. 1: Loschmidt echo for a perturbation in J (solid line) 
with 5.7/ J = 0.05, in U (dashed line) with SU/U = 0.2, and 
for a gravity potential (dash-dotted line) mdg/U = 0.1. For 
comparison, the thick lines show the and dependence. 
All the curves are for a 7 site system with unit density, and 
J/U = 1. Similar behavior is observed for a wide variety 
of initial states. In the inset, schematic of the time reversal 
sequence. The initially prepared state tpo is evolved with the 
optical lattice Hamiltonian for a time t, where the sign of the 
BHM Hamiltonian is changed by a phase imprinting gate and 
a change in the external magnetic field. After another time t, 
the probability of finding the system in tpo is measured, Eq. 
(0. If f{t) = 1, the time reversal was perfect. 



changes drastically the scattering process. In the sim- 
plest cases, the scattering length as a function of the 

magnetic field B is as{B) = ai,g ^1 — Jt%^ ^ , where ai,g 
is the scattering length of atoms in absence of the quasi- 
bound state, Bq is the resonance position (related to the 
energy of the bound state) and AB is the width of the 
resonance. Therefore, to flip the sign of U one needs to 
rapidly change the external magnetic field so that the 
Feshbach resonance is crossed and —as- 

The hopping amplitude J is proportional to the overlap 
between the ground states of neighboring sites 0. Thus, 
in principle, it cannot change sign. However, like the tt 
pulse in the spin echo example above ^„Qi] . one can apply 
an operation that changes the effective sign of the Hamil- 
tonian acting on the wave function. The spectrum of the 
hopping term of the BHM (i.e for U = 0) in each dimen- 
sion is E{k) = — 2Jcos(fcd), where d is the lattice site 
spacing, k — mi/d{N -\- 1) with n = 1..N, and N is the 
number of sites in the lattice. Clearly, the sign of E{k) is 
changed by boosting all k states by momentum Tr/d, i.e. 
E{k + Tr/d) = —E{k). The p momentum translation op- 
erator is diagonal in real space, Pp = exp{—ipd^j jfij): 
this is equivalent to the evolution operator under a lin- 
ear potential F^^jhj for a time r = pd/F. Thus, the 
sign of J is effectively changed by applying a pulsed lin- 
ear potential such that Ft = tt. Alternatively, one can 
understand the change of sign by looking at the effect 
of the linear potential on the Hamiltonian in the Heisen- 
berg picture. Indeed, one can see that aj e^^^aj. 
Thus, the interaction term remains unchanged, and the 
hopping term a\ai+i — > e'-^''a\ai+i\r=-n/F = -ajoi+i- 
In other terms, the linear potential creates a phase im- 
printing on the lattice such that neighboring sites have 
a Ft phase difference. Such a phase imprinting has al- 
ready been demonstrated using a rapid displacement of 
the quadratic trapping potential In general, a pulsed 
laser masked with a linear gradient can be used. The 
linear phase imprint can be done either in one, two, or 
three dimensions. 

The complete sequence for the Loschmidt echo in the 
optical lattice is schematized in the inset of Fig. 1. The 
magnetic field ramp and the tt phase imprinting need to 
be done much faster than the characteristic time of the 
dynamics of the BHM: typically, h/ J ~ 1 ms and h/U ~ 
100 /iS jsj. Notice that it is not necessary to measure 
the full overlap between unknown states. Instead, only 
the probability of finding the system in the initial state 
is needed. I will discuss this issue in detail below, after 
considering the effects of experimental errors in the LE. 

I will classify the sources of decay of the LE in three 
categories: A first category consists of what I call nat- 
ural sources: external degrees of freedom that are not 
described by the simulated Hamiltonian. In this cate- 
gory 1 include e.g. coupling to environments, or trun- 
cated terms of the real Hamiltonian: second neighbor 
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hopping, excitations between different bands, and so on. 
Natural decay sources give a bound on tfie fidelity of 
the quantum simulation. A second category consists of 
artificial sources: external fields or potentials purposely 
placed in the experiment. Since these are not (in princi- 
ple) reversed by the LE, one could measure their strength 
through the decay of the echo. Below I will give an exam- 
ple of how to measure the gravity potential with the op- 
tical lattice. The third category consists of experimental 
errors in the implementation of the LE, which limits the 
observation of decay due to natural and artificial sources. 

The precision oi U ^ —U is limited mostly by the 
homogeneity of the magnetic field in the trapping region. 
Variances across the sample on the order of 10 mG or less 
are possible j^- For comparison, Feshbach resonances 
can be found at fields Bq ~ lOG — lOOOG, with widths 
0{AB) ~ IG. It might be useful to work near a broad 
Feshbach resonance in a region with low sensitivity to 
B. Assuming a homogeneous error SU, and an initial 
Fock state, for short times t <C h/vJSU a perturbative 
expansion of the fidelity f{t), Eq. gives /(t) ~ 1 — 
J^SU^f^. The fourth power contrasts with the typical 
quadratic decay of perturbation theory 0]. It appears 
because the interaction term of the BHM is diagonal in 
the Fock basis, and the second order terms cancel out. 
This expression is valid up to times t oc h/ J, at which the 
second order terms become relevant and a decay f{t) ~ 
exp{—(35U^t^) is observed numerically with /? = 0{U / J) 
, see Fig. 1. 

Perturbation theory predicts that errors in the tt phase 
gate, represented in the BHM by J ^ — J + dJ, cause 
a decay f{t) ~ 1 — SJ'^t'^ for times t <^ h/SJ (see Fig. 
1). Experimentally, the accuracy of the linear phase im- 
printing is limited by fluctuations in the laser intensities, 
the precision of the imprinting mask and the lenses used 
to scale it down to the size of the atomic cloud. Phase 
imprinting resolution of up to ~ 5/Ltm has been achieved 
(compared to lattice spacings ~ 0.5^m), but this is 
for a sharp phase step. Linear smooth grad ients and lat- 
tices with controllable intersite spacing [2l| can make the 
phase imprinting much more precise. Depending on the 
strength of natural or artificial decay sources of the LE, 
precision up to a few percent of J could still be tolerable 
(Fig. 1). The duration r of the phase imprinting pulse 
has to be short enough so that there is no dynamics inside 
each site of the lattice. By expanding a Wannier function 
of width d/iT in the eigenstates of the lineal potential, 
this condition can be cast as t < 2m(P /tt^H. For typical 
values T < 1/xs, within experimental reach (Sl- 

Perhaps the most challenging step of the LE exper- 
iment proposed here is the preparation/detection of a 
particular many-body state. Some states, like the ground 
state of the superfluid phase, or selected Fock states, can 
be faithfully prepared experimentally [Hllll- However, 
measuring the probability of finding the system in one of 
these particular states might not be simple. This could 
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FIG. 2: (left vertical axis) Derivative of the initial decay rate 
a of the LE for two Hamiltonians with J — 5J/2 and J + 5J/2, 
as a function of J/U and fixed SJ = 0.05U. The plots are for 
1-D systems with A'^ = 4, 6 and 8 sites (triangles, squares and 
circles respectively) and one particle per site. The curves are 
normalized to the J = value, qIj^q ~ ~ 1)5 J^. The 
rates are obtained by a Gaussian fit to f{t) = exp{—at^). A 
singularity in the derivative of a develops near the critical 
point of the transition. For comparison, the gap between the 
ground and the first excited state for A'^ = 8 is shown (thick 
line, right axis). In the thermodynamic limit, the transition 
point moves to J = 0.52U [sot - In the inset, the normalized 
decay rates a for the same parameters as the main plot. 



be done for Fock states if the occupancy of each site of 
the lattice can be individually measured: First, the hop- 
ping dynamics is quenched with a sudden increase in the 
optical lattice potential. Then, the number expectation 
in each site is measured, collapsing the state in the Fock 
basis. The process is done many times, with the prob- 
ability of finding each Fock state being its relative fre- 
quency. However, spatial periods of optical lattices are 
shorter than optical resolutions, thus individual site ad- 
dressability is still far from experimental reach. Clever 
techniques 0| or long wavelength lattices could be a so- 
lution to this problem. Also, single site addressability 
might not be strictly necessary. Other schemes with in- 
direct approaches could at least give bounds to the LE. 
For instance, using microwave spectroscopy and density 
dependent transition frequency shifts, the group of Ket- 
terle recently imaged sites in the lattice with different 
occupation numbers Instead of measuring the full 

LE, Eq. 1^, one could prepare and measure a state in a 
subspace of the full Hilbert space, i.e. by removing all 
sites with occupancy 2. The resulting measurement gives 
a bound on the LE determined by the relative size of the 
subspace considered. 

If the BHM is faithfully simulated, the LE can have 
sensing applications. For example, one can measure the 
strength of external potentials {artificial decay sources) 
whose interaction with the atoms does not change sign 
with the LE procedure. The decay rate of the LE mea- 
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sures the strength of the potential p]. Consider a ver- 
tical one dimensional optical lattice, where the atoms 
are affected by an extra linear term given by gravity 
V = dmg'^j juj, with m the mass of the bosons and 
g the gravity acceleration. Time dependent perturba- 
tion calculation gives f{t) ~ 1 — {2gdmJ)'^t'^, supported 
by numerical results with different initial states (see Fig. 
1). The fourth power decay means that performing ac- 
curate measurements of g with the LE could be difficult. 
However, the LE in the BHM might be more sensitive to 
other artificial potentials, and thus be useful for metrol- 
ogy. As a side note, the spectrum of the BHM presents 
many avoided crossings as a function of g, and a regime 
with a Wigner-Dyson statistics of energy spacings. In 
general this is a signature of quantum chaos 1281. The 
LE has been intensely studied in this subject [l3j and 
could be a powerful tool to investigate this problem. 

The LE in the optical lattice can also be used to mea- 
sure the critical point of the superfluid-insulator transi- 
tion in the BHM. The strategy is to use the algorithm 
for detecting quantum critical points with a one qubit 
quantum computer proposed in Ref. However, the 
present approach is simpler to implement because it does 
away with the need of the qubit: the LE is implemented 
directly on the system. The measurement of the critical 
point of the transition would be as follows: First, prepare 
and evolve the ground state of the BHM Hamiltonian 
(this could be relaxed to other states) with a given set 
of parameters U and J. Second, create a "slightly" im- 
perfect time reversal by a known small amount 6, e.g. in 
J. Finally, measure the echo as a function of J, keep- 
ing S fixed. In Fig. 2 the decay rate of the LE for 
short times and its derivative are shown as a function 
of J/U for a fixed perturbation SJ. By increasing the 
size of the system, the derivative of the short time decay 
rate develops a singularity at the critical point, consis- 
tent with the results found in Ref. 17 for a different class 
of systems. In the thermodynamic limit, the 1-D BHM- 
superfluid phase is a critical phase (Kosterlitz-Thouless 
transition). Therefore, the BHM is an interesting sys- 
tem to understand what more information can be ob- 
tained about quantum phase transitions by studying the 
LE Is^. 

In summary, I described how to reverse the time evolu- 
tion of the Bose-Hubbard model that describes ultracold 
bosonic atoms in an optical lattice. Although its realiza- 
tion presents some challenges, it is within reach of current 
or near-future experimental technology. I showed that 
the time reversal (Loschmidt echo) has many interesting 
applications, such as measuring the fidelity of the quan- 
tum simulation of the BHM, measuring the strength of 
external potentials, and even finding the critical point of 
the BHM quantum phase transition. It will be interest- 
ing to see if the LE can be further developed as a sensing 
or a quantum information tool. 
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